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Objective 


Primary: Determine whether RTLS abort at 
T seconds along launch trajectory is 
possible using optimal control theory 


Secondary: Assess effects of bank angle 
constraint, lift coefficient constraint, free 
and fixed final boundary conditions, etc. 


HL-20 PLS BASELINE DESIGN 



29.5 ft r— 


Weight, lb 

Dry (with 22% margin) 19,777 

Landed 22,057 

On-Orbit 26,186 

Launch Escape 8,420 

System/ Adapter 


Removeable panels 



Subsystem 
attachment areas 


Shuttle-derived tiles 


Aluminum 

pressure 

vessel 


L Fighter 
technology 
landing gear 


Gross Launch on NLS 34,607 


HL-20 PLS CURRENT TECHNOLOGY DESIGN 
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THE PERSONNEL LAUNCH SYSTEM (PLS) 


Complementary System to Space Shuttle 

• Space Station crew transfer 
Alternate access to/from space for people/priority cargo 


Space Station Reference Mission 


• Transfer and return up to 
8 Space Station personnel 
and/or priority cargo 

• 72-hour mission duration 

• 1,100 ft/sec on-orbit 
propulsive capability 

• Placed in orbit by existing 
or future booster system 

• Kennedy Space Center 
launch/landing site 

• Alternate landing site 
capability 



Talay 34 


HL-20 Aborts 


VAB Analysis 

• On the pad 

• 0-20 sec 

• 20-65 sec 

• 65-403 sec 

• 403-478 sec 

• 478+ sec 


Return to Launchsite (Shuttle 
landing facility) 

RTLS (Skid strip) 

Ocean landing by parachute 

Transatlantic abort landing 

Abort to orbit 
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Vehicle Aerodynamic Model 


• Aerodynamic data from Jackson and Cruz 


• At each angle of attack and Mach number, find 
$E> 5 L, 5 uthat trim vehicle and minimize drag; 
calculate Cl and Cd here 


• For each Mach number, determine coefficients 
for Cd expression 

C D = C Do (M) + Cd, (M) C l + C Dj (M) C l 


Optimal Control Theory 

• Cost min J = O[X(t 0 ), X(t f )] 

U 

• Plant * = f(%, 0) 

• Constraints: Control State 

g(x, a) = o m) = o 

h(X, 0) < 0 6 (%) < 0 

• Boundary conditions \|/[5t(to), X(tf)] = 0 
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Cost Function and Plant 


• J = -h(tf) (final altitude) 

• States % = [h x y V y \j/] T 

- x,y Cartesian system, x east, y north, origin 

at point runway centerline extended 

- \| f 0 for easterly flight, increases CCW 

• Controls 0 = [C L cy] T 

- a negative for right bank 

• Equations of motion: flat earth, non-thrusting, 

aerospace vehicle 


Control/State Constraints 


• Bank angle can be constrained (40 deg. nominal) 


'^max — C — °max 

• Lift coefficient is constrained between upper 

and lower trim limits (function of Mach) 

C Lml „(M)<C u <C Lm JM) 

• Normal and axial load factor constraints (3 g 

units nominal) 
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Lift Coefficient Constraint 



Mach number 


Initial Conditions 


« Initial conditions for abort at T seconds are 
conditions at T along ascent trajectory 
followed by primary solid rocket motor (srm) 
burn, followed by sustainer srm burn 


Example: Initial conditions for abort at T=30 
h(tO) = 32882 ft V(tO) = 1565 ft/sec 

x(tO) = ■ 7409 ft 7(t0) = 79.7 deg 

y(t0) = 45357 ft \|/(t0) = - 2.0 deg 
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Final Conditions 


28.5 p 

For all cases: 
x(tf) = 0.0 ft 

28 4 - 

y(tf) = 0.0 ft 

Latitude, 

V(tf) = 521 .Oft/sec 1,00 
y(t f ) = -1 9.0 deg 28 3 - 

\|/(tf) = -220.7 deg 

28.2 L 

279.1 



Longitude, deg 


Solution Method 


Trajectory Optimization by Differential 
Inclusion (TODI) 

- eliminates controls from problem by 
constraining state rates 

- leads to nonlinear programming problem 
where parameters are state values at user 
defined nodes (NPSOL) 
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Altitude (ft) 


Solution for 30 Second Abort Case 

Altitude vs. Time 

* 10 ’ 



Solution for 30 Second Abort Case 

Groundtrack 



Solution for 30 Second Abort Case 


Velocity vs. Time 
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o 

o 

O 800 
> 



Time (sec) 


Solution for 30 Second Abort Case 

Flight Path Angle vs. Time 


LL -20 



Time (sec) 


Solution for 30 Second Abort Case 

Heading vs. Time 



-250 H 


0 20 40 60 80 100 120 140 

Time (sec) 


Solution for 30 Second Abort Case 

Lift Coefficient vs. Time 




Solution for 30 Second Case 

Bank Angle vs. Time 



0 20 40 60 80 100 120 140 


Time (sec) 


Comparison of T=30,40,50 Sec. Aborts 

Altitude vs. Time 

* 10 * 
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Comparison of T=30,40,50 Sec. Aborts 

Groundtrack 



Comparison of T=30,40,50 Sec. Aborts 

Velocity vs. Time 






Comparison of T=30,40,50 Sec. Aborts 

Flight Path Angle vs. Time 



0 20 40 60 80 100 120 140 160 


Time (sec) 


Comparison of T=30,40,50 Sec. Aborts 

Heading vs. Time 



0 20 40 60 80 100 120 140 160 


Time (sec) 
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Comparison of T=30,40,50 Sec. Aborts 

Lift Coefficient vs. Time 



Comparison of T=30,40,50 Sec. Aborts 

Bank Angle vs. Time 
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Concluding Remarks 


• When final V is fixed, maximizing final h is 

nearly same as maximizing final energy 
==>calculation of minimum energy trajectories 

• Choice of cost function for abort (and reentry) 

not obvious 

• Future work: 

- Single Stage Vehicle (?) 

- experiment to assess "power” of TODI 
approach compared to traditional shooting 
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Range Optimal Atmospheric Flight Vehicle Trajectories 
in Presence of a Dynamic Pressure Limit 


BY 


Hans Seywald 

Analytical Mechanics Associates Inc. (AMA) 
Spacecraft Control Branch, NASA LaRC 
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Cost Function: 

J [ u ] = - X (tf) 


State Equations: 

E * (tlT - D) ^ 
h = v siny 
y = f (n - cosy) 
x = vcosy 


Initial Conditions: 

E(0) = 38, 029.[m] 
h(0) = 12, 119.[m] 
y(0) = 0.° 
x(0) = 0.[m] 

Final Conditions: 

E(t f ) = 9000.[m] 

Kt f ) — 942. [m] 
y(t f ) = - 11.5° 
x(t f ) be maximized 


Final time: 


Control constraint: 
0 :s rj s 1 

l/ll S /tmax 

State constraint: 

V < VmaxOl) 


t f = 60[sec] 




• Validate optimality of solutions obtained with other methods 

• Use optimal solutions to develop guidance laws based on 
neighboring optimal control 


• Decide on choice of discretization ( e.g. finite elements) 
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tJ**i,_->****> 
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1) V - Vjmx <0* <5 = 1, 

2) v — Vqux = 0, <5 = 1, 


3) v - v n 


0, d m 1, 


M.o 

5/1 

n > 0 from gj(v — v^) — 0 
n < 0 yhwpi ^(v — v »»i) ” 0 


4) v - v^n “ 0, 6 from ^(v - Vgux} = 0, n singular 

5) v - Vgu, = 0, d from ~l?(y ~ Von) = 0, n = n 


6) v — v m < 0, 


<5=0, n = n. 


1 2 different possible control logics are obtained 

6 cases with v^-limit not active 
1 first-order singular case with v max -limit not active 
6 cases with active vw -limit 
1 first-order singular case with v^-limit active 
1 second order singular case with v m , T -limit active 


• To perform higher order optimality tests the Generalized 
Legendre— Clebsch condition has been extended to the case of singular 
control in presence of state/control constraints 
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RE MA RKS 


• Switching structure is non-intuitive 


• Dynamic pressure constraint makes problem very ill-conditioned 

(i) standard shooting codes fail 

(ii) developed flexible shooting code 

(iii) trick: start integration at the end of singular control 





STRUCTURE OF SHOOTING CODE 



• Boundary value problem: find x such that F(x)=0 


• User completely determines function F 


• Simple structure allows independent debugging of F(x) 
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650 
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SUMMARY 



• All possible control logics are analyzed 


• Optimal switching structures are identified. 

Solutions involve singular control along state constrained arcs 


• A flexible multipoint shooting code was developed and applied 
successfully 


• TODI was used to perform sanity check and to guess the optimal 
switching structure 
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3-D Air-to-Air Missile Trajectory Shaping Study 


Renjith Kumar, Hans Seywald 

Analytical Mechanics Associates Inc., Hampton, Virginia 


Eugene Cliff, Late Henry Kelley 

Department of Aerospace Engineering, VPI&SU, Blacksburg, Virginia 



Sir Francis Drake and “Manouevre Board ’ 1 


World War II 

- Pure-Pursuit 

- Deviated Pursuit 

- Command to Llne-of-sight 

- Collision course 

- Proportional Navigation 

Singular Perturbation (Reduced-Order Modeling) 
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u(0 = u N (t,) + G,(r,)[X(/) - X n (/ 7 )] t, t $ t f 

+ [G (t/) + u N (t /)] [f - tj\ + G 2 itj)d\p{t) 
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SALIENT CONTRIBUTIONS 


• Identified attainable sets via intricate homotopy procedures. 

• Checked sufficiency conditions for weak local optimality. 

- Derived a new matrix differential equation for conjugate point 
testing. 

• Developed an efficient method of optimal gain evaluation. 

• Developed a composite midcourse guidance strategy (half-pn) which 
saves on-board storage. 
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Constrained Minimization of 
Smooth Functions Using 
A Genetic Algorithm 


Lynda J. Foernsler, SCB 
Dr. Daniel D. Moerder, SCB 
Dr. Bandu N. Pamadi, Vigyan 


LaRC Workshop 
March 18-19 


Purpose 


• Discuss the use of a simple genetic algorithm 
for constrained minimization of differentiable 
functions with differentiable constraints. 

• Assess the performance of this approach 
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Outline 


• Genetic Algorithms (GA) 

• Problem Formulation for GA 

• Numerical Experiment 

• Comparison to Penalty Function Approach 

• Conclusions 

• Future Work 


Genetic Algorithms 


• Nonderivative, nondescent, random search 

procedures for unconstrained functional 
minimization -* — 

• Algorithmic structure is based on notions 
from biology with "survival of the fittest" 
search heuristic 

• Operations performed on successive 
generations of a population represented 
by binary coded strings (DNA-analog) 
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GA Operations 


M = population 



• Initial population, M 0 , is randomly generated 


Constrained Function Minimization 


x* — min c(x) 

x6«" 


fi(x*) = 0 i € E 
/;(**) >0 jei 
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subject to 







Kuhn-Tucker (KT) Conditions 


where 


dC(x, A) 


dx 


= 0 


X* 


A* > 0 
Kfk(x*) = o 

fk = o 
fk> 0 


j£l 
k e Eul 
keE 
kel 


£(x,\) = c(x)~ Kfk(x*) 

keEui 


Problem Formulation For GA 


• Convert the solution of the necessary conditions 
for a constrained minimum into an unconstrained 
function minimization 

• Solve the resulting unconstrained minimization 
problem 

x* = arg min g(x) 

xe X 


where X is the user-specified bounded volume 
over which the GA takes place. 
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Unconstrained Minimization Problem 

Formulation 


g(x,\*) = ^2 \C x Jx, A*)| + Y] \fj(x)\ + y'|min{0, f k (x)}\ 
i = i jeE kei 

• estimate A* by setting C x (x, A) = 0 




f Vi(x) 
Mx) = ^ \Vi(x)\ 

0 


fz = 0 i e E 

fi< 0 i e I 

fi> o ie I 




KT conditions are satisfied by solving the nonsmooth 
equation 

g(x,v(x)) = 0 


Genetic Algorithm Function Minimization 


x* = argminp(x, v{x)) 

x ^ X 

where X is the user-specified bounded volume over which 
the genetic search takes place: 


X — { X . (^i)min ^ X{ ^ ^ 
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GA Function Minimization 


x * = argmin v(x)) 

where X is the user-specified bounded volume over which 
the genetic search takes place: 


X — { X I ^ ^ (•£ t) max j i — !»•••» 


Numerical Experiment (1) 


• Mission: Determine control settings for an energy-state 
approximation of minimum-fuel ascent to orbit for the 
Langley Accelerator 

• Control variables: 


X = 


' a, 
h, 

\ &E, 

< 5 t , 




angle of attack (deg) 
altitude (ft) 

elevon deflection (deg) 
thrust vector angle (deg) 
fuel equivalence ratio 
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• Cost: 


c{x) = 


dE_ 

dm 


subject to 

- vertical acceleration balance equality constraint 

- pitch moment balance equality constraint 

- dynamic pressure inequality constraint 

• Monte Carlo Experiment 

- 1 00 GA runs 

- 600 generations/run 

• Used final generation x values from GA runs as initial 
guesses for Newton-Raphson (NR) method 


Distribution of KT Error 
for Aerospace Plane Model 



KT Error 
5 27 


Number of Generation 


Distribution of Control Settings 
for KT Approach 


elevon deflection (deg) thrust vector angle (deg) thrust vector angle (deg) 



alpha (deg) alpha (deg) elevon deflection (deg) 


KT Error Thresholds 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

KT Error 
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Number of Runs Satisfying KT Error Thresholds 



KT Formulation Results 


• 82 of the NR runs converged 

• 99/100 runs converged within a KT 
error threshold of .9 

• Fewer number of generations/run would 
have sufficed 


Comparison To Penalty Approach 


Penalty function form: 


z pen = argmm 


c O)+ P( x >fk( x )) 

k<E£ul 


Monte Carlo Experiments 

- 100 GA runs 

- 600 generations/run 

- various penalty-weighting combinations 
Initial Guesses for NR method 5 2 9 


Best Penalty Function Histogram 



0 5 10 15 20 25 30 

KT Error 


Best Penalty Function Scatter Plots 


elevon deflection (deg) 



alpha (deg) 


thrust vector angle (deg) 

20—i 

15 — 





0 12 3 

alpha (deg) 


thrust vector angle (deg) 



-10 0 10 20 
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Penalty Function Results 


• 74 of the NR runs converged for the best case 

• fine tuning of penalty-weighting combinations 
is problem specific 


Conclusions (1) 


• Discussed search characteristics and algorithmic 
operations of a simple genetic algorithm. 

• Discussed method of adapting the KT conditions for a 
constrained minimization problem to formulate an 
unconstrained minimization function to be used by a 
genetic algorithm. 

• Demonstrated KT method formulation numerically on 
an aerospace plane model of the Langley Accelerator 
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Conclusions (2) 


• For this study, KT approach provides 
reliable initial guesses for Newton-Raphson 
method 

• Unlike the penalty approach, the KT 
approach 

- minimizes a function whose optimum 
value is known a priori 

- provides a measure of the constrained 
stationarity of the solution 


Future Work 


• Exploit stopping criterion of KT approach 

• Extend GA algorithm to include non-smooth 
cost function and non-smooth constraints 
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